Lattice and elastic constants of titanium-niobium monoborides containing aluminum and vanadium 
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First-principles electronic-structure computes the lattice and elastic constants of single-crystal TiB and NbB 
and changes with Nb, Ti, Al, and V solutes. The data is built into an interpolation formula for lattice and elastic 
constants of the quartenary (TiNbAlV)B with dilute Al and V concentrations. The lattice and elastic constants 
of borides in two Ti alloys containing Nb and Al are predicted from microprobe measurements. 

Keywords: Elastic behavior; Titanium alloys; Transition metals; Borides; Density functional 



C3 



-a 
c 

o 
o 



> 

00 

m 
o 
o 

-I— > 



I 

c 

o 
o 



X 



Introduction. Strengthening titanium alloys with 
titanium-monoboride (TiB) combines the beneficial stiffness 
properties of TiB with the plastic behavior of a titanium 
alloy matrix, increasing the range of applicability for auto- 
motive and aerospace applications. The Ti-TiB composites 
have increased stiffness, high-temperature strength, creep 
performance, fatigue and wear resistancejll 0, Borides 
also influence the morphology of a phase in a/@ Ti alloys, 
producing microstructures similar to recrystallized alloys|4]. 
Despite the importance of TiB, accurate predictions of 
titanium-monoboride elastic constants have only happened 
recently 1 5]. Titanium-boride has the FeB structure: a primi- 
tive orthorhombic structure with space group Pnma |6j. The 
lattice constants measured with x-ray powder diffraction are 
fl = 6.12 + 0.01 A, b = 3.06 + 0.01 A, and c = 4.56 + 0.01 A 
131. The eight-atom unit cell has four equivalent Ti atoms at 
the c Wyckoff positions (0.177, 0.25, 0.123) and four equiva- 
lent B atoms at the c Wyckoff position (0.029,0.25,0.603). 
Direct measurements of the 9 independent single-crystal 
elastic constants are difficult, and instead elastic-constant 
combinations have previously been inferred from Ti-TiB 
composite measurements! 8, 9]. Panda and Chandran have 
filled the gap of experimental measurements with density- 
functional calculations of the elastic constants of pure 
TiB|H. 

In order to make accurate predictions about the lattice and 
elastic constants of monoboride inclusions in titanium al- 
loys, it is necessary to include the effects of off-stochiometric 
chemistries sampled in the monoborides. These titanium al- 
loys often include aluminum, vanadium, and niobium addi- 
tions. While the borides formed in these alloys are expected 
to remain near 50 at.% B composition, potentially any of the 
four metals may occupy the metallic sublattice. Microprobe 
measurements of borides in titanium alloys containing nio- 
bium have found a nearly equal concentration of Ti and Nb 
in the borides lioll . Metallographic investigation of alloys in 
the Ti-V-B system has found comparable V concentration in 
the metal matrix and borides after partitioning! 1 !]■ Given the 
difficulty of elastic property measurements for the small dis- 
persed borides, density-functional theory can be used to pre- 
dict the missing data. Here, a predictive model of lattice and 
elastic constants for (TiNbAlV)B in the dilute Al and V limit 
is constructed using density-functional theory calculations of 



TiB and NbB properties, combined with misfits due to Nb, Ti, 
Al, and V solutes, all at zero temperature. The predictions 
of the model are compared with a quasirandom (TiNb)B al- 
loys to demonstrate the accuracy of the interpolation. Finally, 
the lattice and elastic constants are predicted for borides using 
compositions taken from microprobe measurements in actual 
allovs lfioll . 

Methods. The ab initio calculations are performed with 
vasp flU fllll . a density-functional code using a plane- wave 
basis and the projector augmented-wave (PAW) methodQ, 
with potentials generated by Kresse|15]. The generalized- 
gradient approximation of Perdew and Wang is employed! 16|. 
In order to ensure accurate treatment of the boron poten- 
tial, a plane-wave kinetic-energy cutoff of 319 eV is used. 
The fc-point meshes for the 8-atom pure monoboride cells are 
7 x 14 x 10 and 7 x 5 x 5 for the 48-atom 1x3x2 super- 
cells, with a Methfessel-Paxton smearing of 0.2 eV. The PAW 
potentials for boron and aluminum treat the s and p states as 
valence, while the titanium, vanadium, and niobium potentials 
treat the s, d, and filled-/? states as valence, corresponding to 
a core atomic reference configurations of He for B, Ne for Al, 
Mg for Ti and V, and Ca for Nb. In all cases, the internal atom- 
istic forces are relaxed to less than 5 meV/A, and the stresses 
are relaxed to less than 0.2 kbar to determine lattice constants. 

The changes in lattice and elastic constants of TiB and NbB 
with dilute substitutions of Nb, Ti, Al, and V are extracted 
from a systematic series of relaxed and strained supercells, 
all at OK. First, the TiB, NbB, and (48-atom) quasirandom 
(Tio.7 5 Nbo.25)B, (Tio. 5 Nbo. 5 )B, and (Tio.2 5 Nbo.75)B lattice con- 
stants are determined via conjugate-gradient minimization. 
Special quasirandom structures 1 1711 1811 are periodic supercells 
where the chemical species are chosen to most closely approx- 
imate a random structure; c.f., Figure ^ The quasirandom 
structures were constructed using the program atat 11911 . To 
determine the elastic constants, relaxed cells are subjected to 
different magnitudes of one volumetric strain and six volume- 
conserving strains (c.f., Table [Q; the magnitudes range from 
5 = -0.005 to 6 = +0.005 in steps of 0.001. After the atomic 
degrees of freedom in each strained cell are fully relaxed, the 
stresses are calculated. Each stress cr,- is a linear combination 
of strains ej and elastic constants Cy: <x, = 2 ; Cy-ey; c.f., Ta- 
ble |U The ratio of stress to 6 from positive and negative 6 
values are averaged. A range of strains is used to check that 
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TABLE I: One volumetric and six volume-conserving strain combi- 
nations and the linearized stress responses used to compute elastic 
constants in monoborides. The magnitude of each strain is given by 
a single parameter 8 that can be positive or negative. The linear re- 
sponse of stress to strain is determined by a linear combination of 
elastic constants, and these linear combinations are inverted to deter- 
mine the elastic constants from the electronic-structure calculations. 
The first strain combination is a purely volumetric strain; the remain- 
ing six combinations conserve volume for all 8. The first four strain 
combinations produce stresses dependent on Cn, Cxi, C33, Cn, C23, 
and C13; the final three give C44, C55, and C^. 
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TABLE II: Lattice and elastic misfits for Al and V substitutions in 
TiB and NbB. The misfit for a lattice or elastic constant y for solute 
X is defined as y~ l dy/dc x (c x = 0). The values are extracted from 
48 atom monoboride supercells, corresponding to an impurity con- 
centration of 4.2 at.% The overall effect of Al and V additions is a 
volume shrinkage for TiB and NbB. Al weakens the elastic response 
of TiB and NbB, while V strengthens the elastic response of TiB and 
weakens NbB. 
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the effect of anharmonicity is negligible; the elastic constant 
combinations change by less than 0.1% between S = 0.005 
and 6 = 0.001. The values for C44, C55, and C(,(, are extracted 
directly from corresponding strains; the other six elastic con- 
stants Cn, C22, C33, C12, C23, and Co are extracted from a 
least-squares fit of the 12 data from the first four strains in Ta- 
bleU To compute lattice and elastic misfits, 1x3x2 supercells 
of 24 metal ions and 24 boron ions are constructed; a single 
Al, V, Ti or Nb ion is substituted for one metal ion in the cell. 
Each 48-atom cell is relaxed to extract lattice constants; the 
cells break orthorhombic symmetry slightly, producing mon- 
oclinic strains of less than 10 -3 . These monoclinic strains are 
removed before calculation of elastic constants, which pro- 
ceeds in an identical manner to the TiB and NbB cases, using 
6 = +0.005 for the volumetric strains, and 6 = +0.004 and 
+0.002 for the volume-conserving strains. The misfit for lat- 
tice or elastic constant y of each metal species X=A1, V, and 
Nb is numerically extracted as 
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for TiB, and similarly for NbB with X=A1, V, and Ti. 

The composition of a general monoboride (TiNbAlV)B is 
given in terms of the atomic compositions of the four metal 
elements, cji, CNb, cai, and cy, under the assumption that 
metal ions only occupy the metal sublattice, and boron ions 
the boron sublattice. The alloy concentration is parameterized 
in the dilute Al and V limit, assuming a random metal sublat- 



tice, with equal partitioning of Al and V to Ti and Nb sites, 
[Ti(i-x)(i-(ff+/j))Nb A (i_( ff+y8 ))Al ff V jg J B : 
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where x is the relative addition of Nb to Ti, and a and j3 are 
the additions of Al and V, respectively. The dilute Al and V 
limit corresponds to a <k 1 and /3 <K 1. Lattice and elas- 
tic constants for a given monoboride alloy concentration are 
interpolated from the limiting TiB and NbB values using the 
cubic polynomial in x, 
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where y represents the given lattice or elastic constant, and all 
derivatives are evaluated in the dilute concentration limit from 
Eqn. 0. The interpolation formula builds the quartenary ran- 
dom alloy response from the six dilute binary alloy responses. 
It assumes that the misfits sum linearly, and interpolates prop- 
erties of (Tii_ v Nb A )B as a cubic polynomial in x. 

Results. Figure |2 and Figure show the lattice and elas- 
tic constants for the random (Tii_ v Nb c )B alloy as interpolated 
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FIG. 1: Special quasirandom structures for (Tio.7sNbo.25)B, (Tio5Nbo.s)B, and (Tio.2sNbo.75)B. The 48-atom supercell of Ti (large dark grey), 
Nb (large light grey), and B (small grey) is repeated in the a direction for visualization. The chemical identity in each periodic supercell is such 
that the neighbors for each atom approximates a truly random structure. These periodic supercell arrangements allow for efficient calculation 
of properties of a random alloy at the given compositions. 
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FIG. 2: Lattice constants for (Ti|_. v Nb v )B as a function of atomic concentration x of Nb. The values for pure TiB and NbB are shown as the 
endpoints of the interpolation curves. The solid line shows the cubic interpolation of Eqn. J3J for a = ft = 0; the straight dashed line indicates 
the result of a purely linear interpolation. The lattice misfits y~'dy(TiB)/dcNb and y _1 dy(NbB)/rffTi are shown italicized with the corresponding 
slopes at the endpoints. For comparison, the lattice constants of the quasirandom (TiNb)B alloys are shown as single points for x = 0.25, 0.50, 
and 0.75. The deviation from the predicted values and the quasirandom alloys are due to non-dilute concentrations reducing the strength of 
misfit. 



from the TiB and NbB endpoints (Eqn. (|3}), and compared 
with quasirandom (TiNb)B alloys. The results for pure TiB 
(x = 0) match previous calculations and experiments Q 01- 
The cubic interpolation of Eqn. (0 captures the effect of mis- 
fit for the small x and (1 — x) limits. There was little difference 
in lattice and elastic constants and total energy between the 
quasirandom structures and the ordered structures, indicating 
that ordering will be weak, and has only a small effect on the 
properties of the borides. The quasirandom (Tio.5Nbn.5)B al- 
loy deviates from the prediction from misfits generally by re- 
ducing the effect of misfit, approaching the average of the TiB 
and NbB properties. In general, the cubic interpolation is not 
dramatically different from the linear interpolation, with the 
notable exception of C22- In this case, the quasirandom alloys 
shows significant deviation in both b and C22 from the mis- 
fit prediction, especially for (Tio.5Nb 5 )B; these are expected 
to be related as decreasing lattice constants are connected to 



increasing associated elastic constants. The large deviation 
indicates that experimental measurements of lattice constants 
and elastic for (Tii_ v Nb x )B may provide potential insight into 
the effect of non-dilute concentrations on lattice and elastic 
properties. 

Table HTl shows the Al and V lattice and elastic misfits in 
TiB and NbB. Generally, Al and V solutes decrease the lat- 
tice and elastic constants of TiB and NbB; the exception is V 
solutes in TiB, which have positive elastic misfits. The reduc- 
tion of elastic constants contrasts the effects of Nb and Ti on 
TiB and NbB, which expand the lattice (c.f., Figure 0; this 
suggests that Al and V additions in the borides offer control 
of lattice matching by offsetting Ti and Nb expansion. The 
decrease in elastic response of TiB with Al and NbB with Al 
and V additions generally matches the Nb and Ti response, 
including the positive misfit for C22; this suggests less possi- 
ble control over elastic response in the borides. While Al and 
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FIG. 3: Elastic constants for (Tii_ v Nb v )B as a function of atomic concentration x of Nb. The values for pure TiB and NbB are shown as the 
endpoints of the interpolation curves. The solid line shows the cubic interpolation of Eqn. J3J for a = ft = 0; the straight dashed line indicates 
the result of a purely linear interpolation. The elastic misfits y~'dy (TiB) /dctib and y~'dy (NbB) /dcji are shown italicized with the corresponding 
slopes at the endpoints. For comparison, the elastic constants of the quasirandom (TiNb)B alloys are shown as single points for x = 0.25, 0.50, 
and 0.75. The deviation from the predicted values and the quasirandom alloys are due to non-dilute concentrations reducing the strength of 
misfit. The general trend is a weakened elastic response in the random alloy compared with the linear interpolation of the endpoints. 



V are present in many Ti alloys, the prediction of segregation 
to the borides with alloy chemistry and heat treatment is dif- 
ficult; inferences from lattice constants are possible. For ex- 
ample, measurements of the lattice constants of the borides in 
Ti-6A1-4V (10 at.% Al, 4 at.% V) give a = 6.11A, b = 3.05A, 
c — 4.56A 1 20] ; assuming a sensitivity of +0.01A, this pre- 
dicts a maximum concentration of 8 at.% Al and 3 at.% V in 
the boride. Generally, experimental measurement of mono- 
boride chemistry is crucial data for the prediction of boride 



lattice and elastic constants in real alloys. 

Table lllllcontains the predicted lattice and elastic constants 
for monoborides in two different Ti alloys, using the mono- 
boride chemistry from microprobe measurements! 10]. The 
experimentally measured monoboride chemistry of the al- 
loys both show nearly equal segregation of Ti and Nb to the 
borides, and appears to be nearly independent of the Nb con- 
centration in the alloy. Al shows a small segregation to the 
boride, which validates the treatment of dilute Al in Eqn. Q- 
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TABLE III: Predicted single-crystal lattice and elastic constants for 
experimentally measured monoboride chemistries] 10]. Microprobe 
analysis of monoborides in Ti-15Al-33Nb-5B (at.%) and Ti-22A1- 
26Nb-5B (at.%) alloys is used to interpolate lattice and elastic con- 
stants using Eqn. As can be seen from Figure [3] the elastic con- 
stants near x = 0.5 are generally reduced below the average response 
ofTiBandNbB. 



chemistry [at.%] 
Ti 22.32 
Nb 29.50 
Al 0.429 
V — 
B balance 



chemistry [at.9i 
Ti 22.23 
Nb 31.46 
Al 0.003 
V — 
B balance 



_Ti-15Al-33Nb- 
lattice [A] 
a = 6.1840 
C n =4345 
Ci 2 = 1376 
C 44 = 1906 
param: x = 0. 

. Ti-22Al-26Nb- 
lattice [A] 
a = 6.1835 
C n =4387 
C l2 = 1396 
C 44 = 1907 
param: x = 0. 



5B (at.%) 

and elastic [kbar] constants 
b = 3.1262 c = 4.6833 
C 22 = 4788 C 33 = 4401 
C 23 = 1237 C 13 = 1133 
C 5S = 1726 C 66 = 2227 
569; a = 0.0082; = 
5B (at.%) 



and elastic [kbar] constants 
= 3.1308 c = 4.6853 
Cn = 4808 C 33 = 4461 
C 23 = 1242 Ci 3 = 1136 
C ss = 1756 C 66 = 2249 
586;a = 5.6xl0- 5 ;/? = 



In both cases, the final B concentration is nearly 50 at.%; 
this is expected from the Ti-B phase diagram, which shows 
TiB to be a line compound with a homogeneity range of 49- 
50 at.% ItI EUE^ . Microprobe measurements do not provide 
information regarding the possibility of local short-range Ti- 
Ti or Nb-Nb order in the boride. As Figure[5]shows, the lattice 
constants are expanded relative to the linear interpolation of 
TiB-NbB; Figure [3] also generally shows a reduction in elas- 
tic constants relative to the linear interpolation of TiB-NbB. 
Experimental verification of the predicted lattice and elastic 
constants can elucidate the effect of non-dilute concentrations 
in the boride lattice, through deviations from the random-alloy 
predictions. 

Conclusion. The results illustrate the ability of electronic- 
structure calculations to predict lattice and elastic constants 
for real materials, especially in cases where experimental 
measurements are extremely difficult or time-consuming. The 
calculations for (Tii_ A Nb v )B borides and the inclusion of di- 
lute Al and V solutes provides a database for the prediction of 
lattice and elastic constants in real Ti alloys, which in turn can 
aid in the design of new materials. To model increasingly re- 
alistic material systems requires modern computational mate- 
rials science techniques that investigate the changes in chem- 
istry produced by solid solution, and the effect on material 
properties. 
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